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DISTRIBUTED  SYSTEM  OPTIMAL  CONTROL  AND  PARAMETER  ESTIMATION:  COMPUTATIONAL 
TECHNIQUES  USING  SPLINE  APPROXIMATIONS 


H.  T.  Banks 

Lefschet:  Center  for  Dynamical  Systems,  Division  of  Applied  Mathematics, 
Brown  University,  Providence,  RI,  USA 


Abstract .  Spline-based  computational  procedures  for  parameter  estimation  and 
optimal  control  problems  involving  delay  and  partial  differential  equations 
are  outlined.  Brief  discussions  are  presented  reporting  on  application  of 
of  the  ideas  to  compute  feedback  controls  for  delay  systems,  to  estimate 
variable  coefficients  in  population  dispersal  models,  and  to  estimate  para¬ 
meters  in  higher  order  models  arising  in  elasticity. 

i eywords .  Parameter  estimation;  feedback  controls;  delay  and  distributed 
systems;  computational  methods;  splines;  wind  tunnels;  ecology;  elasticity. 


INTRODUCTION 

In  this  note  we  present  a  brief  summary  of 
some  computational  methods  that  have  proven 
useful  in  optimal  control  (both  open  loop  and 
feedback)  and  parameter  estimation  problems 
for  certain  distributed  systems.  We  first  out¬ 
line  the  basic  ideas  which  are  common  to  our 
approach  whether  we  are  dealing  with  function¬ 
al  differential  equations  (FDE)  or  partial 
differential  equations  (PDE)  of  hyperbolic  or 
parabolic  type.  Roughly  speaking,  in  each 
case  one  views  the  system  under  consideration 
as  an  abstract  system 

C(t)  =artq)z(t)  ♦  F(q,t) 


in  the  examples  discussed  below).  One  then 
approximates  the  system  (1)  (and  the  corres¬ 
ponding  control  or  estimation  problem)  in  the 

subspace  Z  by  the  system 

:N(t)  -j/'(q)zN(t)  *  PSF(q,t) 


where  P^  is  the  canonical  orthogonal  projec¬ 
tion  of  Z  onto  ZN  and  =  P!]'a*'r''' .  This 
results  in  an  approximate  estimation  or  con¬ 
trol  problem  entailing  a  finite  dimensional 
state  space  to  which  standard  computational 
packages  can  be  applied. 


in  an  appropriatelv  chosen  Hilbe’-t  space  Z. 

In  the  event  that  the  operator  2#  depends  on 
parameters  q  to  be  estimated,  one  then  has 
data  or  "observations",  say  c(t^),  and  one 

attempts  to  choose  a  parameter  q  from  an 
admissible  set  Q  so  as  to  yield  a  best  fit 
of  the  model  (1)  to  the  data.  For  optimal 
control  problems,  the  parameters  q  are  pre¬ 
sumed  fixed  and  known  and  F  in  (1)  is  a 
control  input  term,  say  F(t)  =  Bu(t).  Then 
one  has  a  performance  measure  J  depending 
on  i  and  u,  and  an  admissible  control  set 
“Sr  (either  open  loop  or  feedback).  One  seeks 
a  u  in  ‘2'’  that  minimizes  J  subject  to  (1). 

For  either  FDE  or  PDE  systems,  these  problems 
involve  infinite-dimensional  state  systems 
and  hence  computational  schemes  must  be  based 
on  some  type  of  approximation  idea.  The 
approach  we  describe  here  entails  choosing 
a  sequence  of  finite  dimensional  subspaces 

ZN  of  Z  generated  by  basis  elements  con¬ 
sisting  of  splines  (linear,  cubic,  or  quintic 


The  fundamental  convergence  theory  which  can 
be  used  in  either  control  or  parameter  esti¬ 
mation  problems  is  based  on  semigroup  approx¬ 
imation  results.  Briefly  the  ideas  are  as 
follows  (details  differ  depending  upon  whether 
one  is  treating  FDE  or  PDE).  One  first  demon¬ 
strates  that  £/  (q)  satisfies  a  unifo~m  dis¬ 
sipative  inequality  in  Z  (such  as^fqlc.cS 
_<  u>  for  z  €  Dorn  (of{q))i  and  ja^(q) 

(or  some  maximal  dissipative  extension)  gen¬ 
erates  a  Cp- semi  group  T(t;q).  The  approxi¬ 
mating  operators  (q)  are  defined,  as  we 

have  already  indicated,  byjj^  (q)  =  P^af;q)P^, 
and  generate  a  stable  family  of  schemes  such 

that  |exp|U('^(q)t]  |  _<  Mewt  where  M  and 
tv  are  independent  of  q  and  N.  One  uses 
standard  estimates  from  spline  approximation 

theory  to  argue  that  (q):  -*  in 

an  appropriate  sense.  One  then  employs  the 
Trotter-Kato  theorem  (a  functional-analytic 
version  of  the  Lax  Equivalence  theorem: 
stability  plus  consistency  yield  convergence) 

to  establish  that  explja^  N(q)t]  T  ( t ;  q ) 
strongly  in  Z,  and,  moreover,  that 


V*. 


-v*-‘ 
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X 


N 


lt;q)  ;(t;q)  where  ;  and  :  are  solu¬ 
tions  of  (2)  and  (1),  respectively.  These  con¬ 
vergence  results  are  used  in  turn  to  argue 

-Nh 

existence  of  an  estimate  q=lim  q 
for  the  parameter  estimation  problem  involving 
_Ng  .V; 

(1) .  where  {q  }  is  some  subsequence  of  (q), 

_  \ 

q'  a  solution  of  the  estimation  prohelm  for 

(2) .  In  the  case  of  optimal  control  problems, 
one  uses  the  convergence  results  to  establish 


-Ni 
u  k 


to  u,  where 


convergence  of  a  subsequence 

u^  and  u  are  solutions  of  the  control  prob¬ 
lems  for  12)  and  (1),  respectively. 

The  necessary  basic  spline  approximation  theory 
for  FOE  can  be  found  in  (Banks  and  Kappel,  1979), 
whereas  applications  of  the  ideas  to  parameter 
estimation  and  control  problems  are  given  in 
(Banks,  Burns,  and  Cliff,  1981),  (Banks  and 
Daniel,  1981a),  (Daniel,  1982).  Fundamental 
theory  for  PDE  estimation  and  control  prob¬ 
lems  is  developed  in  (Banks  and  Kunisch,  1981) 
and  (Banks,  Crowley,  and  Kunisch,  1981).  Ke 
shall  not  discuss  further  the  theoretical 
aspects  of  our  approach  in  this  presentation; 
rather  we  shall  briefly  outline  several  speci¬ 
fic  problems  in  which  use  of  these  methods  has 
been  quite  fruitful.  The  first  application 
involves  optimal  feedback  control  for  a  delay 
system  problem,  while  the  second  and  third 
deal  with  parameter  estimation  problems  for 
PDE . 


FEEDBACK  CONTROL  FOR  DELAY  SYSTEMS 

The  problem  of  constructing  feedhack  controls 
for  hereditary  or  delay  systems  is  not  new  and 
there  is  a  rather  large  literature  which  we 
shall  not  discuss  here.  Our  own  renewed  interest 
in  this  problem  was  motivated  by  problems 
arising  in  the  design  of  controllers  for  a 
liquid  nitrogen  wind  tunnel  (the  National 
Transonic  Facility  or  NTF)  currently  under 
construction  by  NASA  at  its  Langley  Research 
Center  in  Hampton,  VA.  Kith  this  wind  tunnel 
it  is  expected  that  researchers  will  be  able 
to  achieve  an  order  of  magnitude  increase  in 
the  Reynolds  number  over  that  in  existing 
tunnels  while  maintaining  reasonable  levels 
of  dynamic  pressure.  Test  chamber  temperatures 
(the  Reynolds  number  is  roughly  inversely 
proportional  to  temperature)  will  be  maintained 
at  cryogenic  levels  by  injection  of  liquid 
nitrogen  as  a  coolant  into  the  airstream  near 
the  fan  section  of  the  tunnel.  In  addition  to 
a  gaseous  nitrogen  vent  to  help  control  pressure, 
motor  driven  fans  will  be  used  as  the  primary 
regulator  of  Mach  number.  Fine  control  of 
Mach  number  will  be  effected  through  changes 
in  inlet  guide  vanes  in  the  fan  section. 
Schematically,  the  tunnel  can  be  depicted  as 
in  Fig.  1 . 
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Fig.  1. 

The  basic  physical  model  relating  states  such 
as  Reynolds  number,  jressure,  and  Mach  number 
to  controls  such  as  LN,  input,  GX,  bleed, 

and  fan  operation  .nvolves  a  formidable  set  oi 
PDE  (the  Navier-Stokes  theory)  to  describe 
fluid  flow  in  the  tunnel  and  test  chamber. 

This  model  has,  not  surprisingly,  proved  to  be 
very  unwieldy  from  a  computational  viewpoint 
and  is  difficult,  if  not  impossible,  to  use 
directly  in  the  design  of  sophisticated  control 
laws.  (Both  open  loop  and  feedback  controllers 
are  needed  for  efficient  operation  of  the 
tunnel  -  and  this  is  a  desirable  goal  since 
cost  estimates  for  liquid  nitrogen  alone  arc 

$b.F>  *  1CI6  per  year  of  operation.  '  In  ajd.- 
tion  to  the  design  of  both  open  loop  and 
closed  loop  controllers,  parameter  estimation 
techniques  will  be  useful  once  data  from  tin 
completed  tunnel  is  available  (current  nice- 
tigations  invoice  use  of  data  from,  a  1  5 
meter  scale  model  of  the  tunnel  I . 

In  view  of  the  schematic  in  fig-  1,  it  i-  not 
surprising  that  engineers  le.g.  see  (Armstrong 
and  Tripp,  19811  and  (Gumas,  19"91)  have  pro¬ 
posed  design  of  control  laws  for  subsvs terns 
modeled  by  lumpeJ  parameter  models  (the 
variables  represent  values  of  states  and  con¬ 
trollers  at  various  discrete  locations  in  the 
tunnel  and  test  chamber i  with  transport  do  la's 
to  account  for  flow  times  in  sections  of  the 
tunnrl  A  specific  example  is  the  model 
(Armstrong  and  Tripp,  1981'  for  the  Mach  no. 
control  loop  in  which  variations  in  the  Mach 
no.  (in  the  test  chamber <  arc, to  first  Older, 
controlled  h>  variations  in  thr  inlet  guide 
vane  angle  setting  (in  the  fan  section)  i.c. 
,‘M'ti  —  A- i  t  ri  where  r  represents  a  trans¬ 
port  time  from  the  fan  section  to  thi  tost 
section  Mert  precisely,  the  proposed  equation 


1 
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relating  the  variation  4m  (from  steady  state 
operating  values)  in  Mach  no.  to  the  variation 
66  in  guide  vane  angle  is 

1 6*M ( t )  ♦  6M(t)  =  k  6b  (t-r) 

while  the  equation  relating  the  guide  vane 
angle  variation  to  that  66^  of  an  actuator 

is  ..  . 

56  (t)  *  24*66  (t)  +u/(56(t)  =  </  66A(t). 

Rewriting  the  system  in  vector  notation,  one 
thus  finds  that  the  Mach  no.  control  loop  in¬ 
volves  a  regulator  problem  for  the  equation 

x(t)  =  ApX(t)  ♦  AjX(t-r)  +  B0u(t)  (3) 


where  x  =  (6M,  66,6b).  u  =  66  Here  the 

A 

control  is  the  guide  vane  angle  actuator  in¬ 
put.  A  similar  4-vector  system  problem  can 
be  formulated  in  the  case  where  one  treats 
the  actuator  rate  66^  as  the  control  - 

see  (Armstrong  and  Tripp,  1981),  (Daniel,  1982). 

Problems  such  as  that  just  outlined  led  us  to 
consider  the  spline  techniques  of  (Banks  and 
kappel,  1979)  for  computation  of  feedback 
controls  in  regulator  problems  governed  by 
n-vector  systems 

x(t)  =  L(xt)  *  Bu(t)  (4) 

(X(0),  x0)  =  CH0),<P) 

where  x.  denotes  the  function  s  -*  x(t*s), 

t  v  • 

-r  <  s  <  0,  and  L(x.)  =  I  A.x(t-r.)  * 

[0  ~  1  J*0  J  J 

A(s)x(t>s)ds  with  0  »  rQ  <  Tj  <  ...  <  ry « 

r.  The  cost  functional  is  the  usual  integral 
quadratic  payoff 


J(z0,u)  =  J ( (x (0) , xQ) ,u)  - 

r 

j  x(t)PQx(t)  *  u(t)Ru(t) 

where  t  x(t)  is  the  solution  of  (4).  As 
is  well-known  (see  the  summary  and  references 
to  previous  literature  in  (Gibson,  1980))  the 
appropriate  state  feedback  control  is  given 
in  terms  of  a  functional 


u(xt)  =  -R'1BT(K0x(t) 


,0 

K, (s)x(t*s)ds]  (3) 
J-r  1 


where  the  gains  kQ,  kj  satisfy  certain  Riccati 

type  equations.  A  detailed  explanation  of  use 
of  the  spline-based  methods  for  computations 
in  these  problems  is  given  in  (Banks  and  Rosen, 
1982);  we  only  outline  the  procedures  here  and 
discuss  our  numerical  findings  for  the  NTF 
example. 


with  J  =  J(z0>u)  =<’D:0,z0'>  where  the 

bounded  linear  operator  11:2  -►  2  is  the  solu¬ 
tion  of  the  Riccati  algebraic  equation  (RAI  I: 

j*'*  n+  ij/-  z0R~lSi  2  *  CJ>  =  n  . 

Here  Si  and  3>  are  operators  on  2  given  b 
Si{r<  ,$)  =  (61,0)  and  S»(b,$  j  =  (p  '  ,  0  . 
Carrying  out  the  approximation  steps  outlined 
above  (we  used  the  first  order  spline  sub¬ 
spaces  explained  in  detail  in  54,  p .  idly- Si 2 
of  (Banks  and  kappel,  1979)  for  the  calcula¬ 
tions  on  the  NTF  example  and  a  number  of  otlur 
test  examples  reported  on  in  (Banks  and  kc-ei;, 
1982)) one  obtains  an  approximating  RAE  cor¬ 
responding  to  the  associated  regulator  problem 
N 

for  (2)  in  2  .  This  is  given  by 

s*fN‘n .  r  r.  at N-  7  =  o 

N  N  N  N 

where  J*' N  =  P^PN,  Si*  =  P KSi ,  <&*  =  PN^rN. 

This  equation  has  a  matrix  representation  and 
hence  standard  techniques  -  we  used  the  ORACLS 
package  of  Armstrong  (1980)  in  our  calculations 
-  can  be  employed  to  calculate  and  the 

resulting  approximate  optimal  feedback  control 
uN(t)  =  -R'1BT[D°Px(t)  *  j°n‘°(s)x(fs)ds]  (7) 

where  D^°,  are  components  of  2^  which 

approximate  the  gains  Kp,  kj  of  (5). 

We  investigated  the  Mach  no.  control  loop 
problem  of  (Armstrong  and  Tripp,  1981)  de- 
scribedabove  with  the  spline-based  methods 
and  with  the  so-called  "averaging"  approxima¬ 
tion  method  (a  "zero-order"  spline  type  scheme- 
which  can  also  be  developed  in  the  context  of 
Walsh  function  approximations  -  that  is  not 

of  the  form  =  P '&?'  but  nonetheless  can 

be  used  to  compute  gains  in  delay  system  prob¬ 
lems  -  for  details  see  (Gibson,  19801).  For 
this  example  the  matrices  in  (3)  have  the  font, 


with  the  parameter  values  given  hv  X  -  1.9»>4 
** 

sec,  r  =  .33  sec,  =  3b,  4  -  .8,  k  =  -.011'. 


Briefly  then,  or--  reformulates  the  system  (4) 
as  an  abstract  system  (1)  in  the  Hilbert  space 

2  *  Rn  x  l2(-r.0;Rn)  with  z(t)  =  (x(t),xt). 

The  optimal  feedback  control  is  then  given  by 
(see  the  summary  in  §4  of  (Gibson,  1980)) 

u(t)  =  -R'1#  II  z(t)  (6) 


Calculations  were  carried  out  for  the  problem 
of  driving  6M  from  -.1  to  0.0  (corresponding 
to  M  varying  from  .8  to  .9)  and  66  from 
8.55  to  0.0  (corresponding  to  the  guide  vane 
angle  varying  from  10.48°  to  a  steady  state 
of  1.93°).  Excellent  results  were  obtained 
even  for  low  values  of  the  approximation  index 
N  (N  =  2,4,8).  The  corresponding  optimal  con¬ 
trols  (7)  appear  to  converge  rapidly  to  an 
optimal  control  of  the  form  (b)  (of  course,  we 


.  ^cv» 


do  not  know  n  for  this  example)  and  when  used 
as  feedback  in  (3),  yield  trajectories  as 
graphed  in  Fig.  2.  (In  Fig.  2  we  compare 
these  with  similar  trajectories  obtained  in 
(Armstrong  and  Tripp,  1981)  where  finite 
difference  techniques  -  assuming  piecewise 
constant  controls  -  were  used  to  discretize 
the  delay  system  before  applying  standard 
regulator  theory  for  discrete  systems  to  cal¬ 
culate  the  feedback  controllers.) 


TIME  (SEC) 

Fig.  2. 


The  spline  method  was  computationally  more 
stable  and  in  general  more  accurate  than  the 
averaging  method.  Doth  methods  appear  to 
offer  higher  order  approximations  than  the 
finite  difference  techniques  of  (Armstrong 
and  Tripp,  1981).  All  three  methods  appear 
to  be  adequate  for  the  simple  NTF  example 
we  investigated. 


for  detailed  remarks  on  the  computations  along 
with  some  comments  on  theoretical  aspects  (see 
also  (Gibson,  1980)  and  (kunisch,  1980/  for 
reiazed  discussions). 

PARAMETER  ESTIMATION  IN  POPULAl ION 
DISPERSAL  MODELS 


An  important  problem  to  population  ecologist- 
involves  modeling  of  dispersal  or  movement  of 
insects  and  other  herbivores  in  vegetation 
patches.  There  is  a  growing  literature 
(Levin,  lQ'M,  1981  I  ,  (Okubo,  1980,,  (Karel  .... 
1982)  on  the  use  of  mathematical  models,  e- - 
pecially  ones  entailing  distributed  systems , 
to  investigate  the  effect  of  various  types  of 
transport  mechanisms  on  overall  population 
movement.  Many  important  ecological  investi¬ 
gations  result  in  the  need  to  solve  an  inverse- 
problem  for  a  diffusion  or  more  general  trans¬ 
port  equation.  Roughly  speaking,  given  data 
describing  changing  population  densities  of 
a  specie  (or  species)  and  some  preliminary 
hypotheses  about  boundary  conditions  and 
transport  functions  (spatially,  temporally  or 
perhaps  even  density  dependent)  in  a  proposed 
model,  one  desires  to  estimate  or  identify 
parameters  (including  the  transport  functions) 
in  the  model  and  quantify  the  success  (or  lack 
thereof)  of  the  model  in  deicribing  the  data. 
For  example,  typical  models  might  involve  the 
general  transport  equation  (Okubo,  1980,  p.9S' 
(in  one  spatial  variable)  for  population  den¬ 
sity  u  (here  we  mention  only  single  specie 
models  but  coupled  equations  for  multiple 
species  models  could  also  be  treated  with  the 
ideas  we  outline)  given  from  mass  balance 
considerations  by 

U  ♦  i  <Vu>  '-81 


where  the  "directed  movement”  or  advective 
(convective)  term  ~  (Vu)  contains  a  spa- 

dX 


tially  varying  "velocity"  term  V  =  V ( x l .  The 
diffusion  term  is,  as  usual,  a  result  of 
assuming  Lick's  first  law  of  diffusion,  while 
f  represents  a  general  birth/death  term.  In 
such  models  it  is  often  important  to  allow  the 
transport  functions  @  and  V  to  vary  spa¬ 
tially,  temporally,  or  even  with  population 
density  (ot  perhaps  some  combination  of  these). 
Other  basic  transport  assumptions  or  hypotheses 
(e.g.,  see  (Okubo,  1980,  p. 84-88),  (Pobrhansky 
and  colleagues,  1979))lead  to  different  models, 
but  in  most  cases  a  very  important  problem 
consists  of  using  field  data  *  estimate  the 
transport  functions  (such  as  2  and  \  I  and 
perhaps  birth/death  parameters  in  f. 


We  also  tested  the  spline  method  (and  compared 
it  with  the  averaging  method)  on  a  number 
of  other  delay  system  regulator  examples  and 
our  findings  are  detailed  in  (Banks  and 
Rosen,  1982).  In  summary  we  found  the  spline 
method  generally  at  least  as  good  as  and  in 
some  cases  superior  to  the  averaging  method 
with  regard  to  computational  stability, 
accuracy,  and  rate  of  convergence.  We  refer 
interested  readers  to  (Banks  and  Rosen,  1982) 


We  have  successfully  applied  the  spline  method? 
outlined  above  in  connection  with  (1)  and  (21 
to  such  problems  (Banks  and  Kareiva,  1982). 

In  addition  to  (Banks  and  Kareiva,  1982),  one 
may  consult  (Banks,  1981),  (Banks,  Crowley, 
and  Kunisch,  1981)  for  the  theory  behind  our 
efforts.  Briefly,  one  rewrites  (8)  in  the 
form  (1)  in  the  Hilbert  space  Z  =  L-,(0,1) 
and  then  uses  the  approximating  equafion 
(2)  -  in  this  case  we  employed  cubic  splines 
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for  the  basis  elements  in  I'  -  with  the  data 
to  estimate  the  unknown  transport  functions. 
For  the  resulting  finite  dimensional  problems 
we  employed  a  standard  1MSL  package  (ZXSSQ) 
for  the  Levenberg-Marquaidt  algorithm  in  our 
parameter  search  for  a  fixed  level  of  approx¬ 
imation  N.  In  the  particular  problems  we 
investigated,  we  hypothesi zed  equation  (8) 
in  which  V  =  V(x),  ^  is  constant,  and  f 
contains  piecewise  linear  (in  u)  terms  with 
spatially  dependent  coef f 1 ci en t s .  We  also 
hypothesi ted  unknown  parameters  in  the  initial 
population  density.  Our  early  efforts  with 
field  data  collected  by  P.  kareiva  (the 
experiments  involved  the  dispersal  of  flea 
beetles  in  cultivated  collard  patches)  re¬ 
vealed  that  models  such  as  (8)  with  a  spa¬ 
tially  dependent  V  yield  significantly 
better  fits  to  the  data  than  do  models  with 
V  vanishing  or  chosen  as  some  nontrivial 
constant.  Our  more  recent  efforts  (detailed 
in  (Banks  and  kareiva,  1982)),  again  using 
the  flea  beetle  data,  involve  the  particular 
equation  (obtained  from  (8)  after  some  trans¬ 
formations  and  assumptions) 


3u 

3T 


=  9 


du 


3x 


for  t  '  0,  0  £  x  <  1 . 
assumed  to  have  the  form 


q3(x)u  ♦  g(t,x) 
The  function  q7 


(9) 


here  will  also  prove  useful  in  these  problems. 
PARAMETER  ESTIMATION  IN  ELASTIC  STRUCTUM.S 


We  turn  finally  to  a  brief  discussion  of  use 
of  cubic  and  quint ic  spline  schemes  for  para¬ 
meter  estimation  problems  arising  in  the  stud, 
of  elastic  and  viscoelastic  bodies.  Our 
interest  in  such  problems  was  motivated  b;. 
discussions  with  NASA  engineers  who  desire 
to  estimate  material  properties  for  large 
space  structures  from  observations  of  the 
motions  of  these  structures.  The  sing  lest 
components  of  these  antennae  and  space  st.i 
tions  are  beam- like  and  made  from  coir.p-e>  1 1 1 
materials  (e.g.  graphite  epoxy*.  Thus  as  a 
basic  problem,  we  (Banks  and  Crowley,  1961 
considered  estimation  in  equations  such  as 
those  arising  in  the  Eul er-Bernoul 1 i  theor. 
for  transverse  vibrations  of  a  thin  elastic 
or  viscoelastic  beam  which  is  possibly  sub¬ 
ject  to  damping.  More  precisely,  the  well- 
known  equations  for  the  transverse  vibrations 
of  a  thin  elastic  beam  (no  damping!  are 


Jt* 


c'X*’ 

3  •>  -  ,  , 

— ^  =  f  i.  q ,  t ,  x ) 

3x~ 


q;W 


5-l<  x<  .5  *  L 
otherwise, 


where  T  <  0,  and  q,  contains  an  appropriate 
death  rate  term  (modest  death  rates  within  the 
vegetation  patch,  high  death  rate  outside  the 
patchl  in  additionto  a  term  involving  q 2,  . 

The  function  g  contains  terms 


arising  from  standard  transformations 
of  (8)  with  nonhomogeneous  boundary  conditions 
to  (9)  with  homogeneous  conditions.  The  cubic 
spline-based  estimation  techniques  have  per¬ 
formed  extremely  well  in  our  efforts  to  esti¬ 
mate  q^,  q,,  q^  as  well  as  the  initial  con¬ 
ditions  from  the  data.  The  methods  were  very 
stable  and  rapidly  convergent,  yielding  sat¬ 
isfactory  estimates  for  low  (N  =  8,161  values 
of  the  approximation  index. 


There  is  strong  evidence  (Dobthansky  and 
colleagues,  1979)  ,(Aikman  and  Hewitt,  19"2) 
of  the  need  in  certain  population  studies  to 
estimate  time  dependent  transport  coefficients. 
Our  cubic  spline  methods  can  be  developed  for 
these  problems  (see  (Banks  and  Daniel,  1981b) 
for  preliminary  theoretical  results)  and  we 
are  currently  pursuing  investigations  along 
these  lines. 


There  are  also  numerous  important  control 
problems  arising  in  the  context  of  ecological 
investigations.  Once  an  adequate  model  is 
developed  (the  parameter  estimation  problem) , 
one  might  wish  to  estimate  (calculate)  the 
optimal  vegetation  density  in  a  patch  in 
order  to  hold  population  levels  in  the  patch 
to  a  minimum,  or  at  least  below  some  given 
level.  We  believe  that  the  methods  discussed 


where  Jt  is  the  bending  moment,  m  is  the 
mass  per  unit  length  and  f  is  the  applied 
load.  Two  types  of  damping  are  included  in 
our  formulation.  The  first  is  simply  viscous 
damping  Tut  while  the  second  is  structural 

damping.  For  a  Voight  solid  Ithe  simplest 
viscoelastic  model)  one  has  the  constitutor 
relationship  c  =  EE.  *  cc  .  Thus  the  stress 
o  is  no  longer  proportional  to  the  strain 
alone  (as  in  Hooke's  law!  but  a  term  propor¬ 
tional  to  the  strain  velocity  is  added.  Ir 
this  case  the  usual  Euler-Bernoul )i  formula¬ 
tion  becomes 


Jt  = 


r 

j 


C  ydA 


3x 


This  results  in 
3 

mu  *  — t  (El 

3.x  “ 


the  equation 

-i 2  *,5 

o  U  ,  r  u 

- V  ♦  Cl  — tc~ 

3x“ 


1  ♦  HI  =  f.  1 


which  can  be  rewritten  in  the  form  ( 1 ’  m  the 

Hilbert  space  I  =  H~  x  L,  .  We  have  developed 
and  tested  estimation  schemes  i.for  estimation 
of  parameters  such  as  El/m,  cl/m,  Vm)  using 
cubic  and  quintic  spline  approximation  suh- 

spaces  2  modified  to  treat  various  impor¬ 
tant  boundary  conditions  (simply  supported, 
cantilevered,  as  well  as  beams  with  applied 
moments  at  one  end) ,  The  methods  proved 
extremely  efficient  as  the  detailed  presenta¬ 
tion  in  (Banks  and  Crowley,  1981)  documents. 

While  the  Euler-Bernoul 1 i  equation  (ID)  is 
applicable  in  many  applications  (especial ly 
large  space  structures),  a  somewhat  more 
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involved  analysis  is  required  in  situations 
where  rotatory  inertia  and  shear  effects  play 
a  significant  role  in  the  dynamics.  This 
theory  is  often  necessary  when  high  frequency 
oscillations  of  the  beam  must  be  considered 
(e.g.  in  aerodynamic  structures).  In  this 
event  the  Timoshenko  formulation  is  more 
appropriate.  This  theory  can  be  embodied  in 
a  single  higher  order  equation  (fourth  order 
in  t  and  x  derivatives)  where  the  boundary 
conditions  for  even  the  simply  supported  beam 
involve  second  order  derivatives  in  both  x 
and  t.  For  our  purposes  it  is  much  more 
desirable  to  treat  a  system  of  lower  order 
equations  with  the  corresponding  boundary 
conditions.  The  equations  modeling  transverse 
vibrations  of  a  homogeneous  isotropic  elastic 
beam,  including  rotatory  inertia  and  shear 
effects,  can  be  written  in  terms  of  the 
transverse  displacement  y  and  the  angle  v 
of  rotation  of  the  beam  cross  section  from 
its  original  vertical  position  as 


in  which  not  only  the  parameter  L  ,n 
u  =  (Eu  I  ,  but  also  parameters  l  k. 

1 1  x  x  *  i  : 

in  elastic  (u  (t,(l)  ♦  kjU(t,0j»0»  and 

absorbing  (u^(t,l)  ♦  k-,uj.(t,li  -  c  i  boundar* 

conditions  must  be  identified  from  data. 

A  general  theory  plus  numerical  results  ob¬ 
tained  when  applying  our  approximation 
methods  to  nonlinear  hyperbolic  and  parabolic 
PDt  can  be  found  in  (Banks  and  Kunisch,  IbS; 
(Banks,  Crowley,  and  Kunisch,  JP81).  Other 
areas  of  application  in  which  we  have  a  Ire. id. 
used  or  are  currently  using  these  spline 
based  methods  include  estimation  problems  for 
transport  of  labelled  substances  in  brain 
tissue,  determination  of  static  antenna  con 
figuration  and  shape,  and  estimation  of  poro¬ 
sity  and  permeability  in  porous  media. 
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>tt  =  a‘(>x*  '  yx> 


V-  =  b"V  *  c"(v  -VI 
tt  XX  X 


(11) 


with  the  boundary  conditions  for,  say,  a 
fixed  end  beam  given  by  y ( t , 0 )  =  y(t,l)  =  0, 

v(t,0)  =  V(t,l)  =  0.  Here  a2  «  k 1  AG/m, 

*»  ■>  -» 

b*  =  El/m,  c“  =  Aa'/l  with  A  =  cross  sec¬ 
tional  area,  E  =  Young's  modulus,  C  -  shear 
modulus,  I  =  moment  of  inertia,  and  k’  = 
shear  coefficient. 
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Equation  (11)  can  be  rewritten  in  the  form 
(1)  in  the  Hilbert  space  3  *  H*  *  L.,  x  H*  x 

L,  and  then  cubic  spline  schemes  can  be 
applied  (the  approximating  equations  again 
have  the  form  (3))  to  estimate  parameters 
such  as  a,  b,  and  c.  We  did  this  (Banks 
and  Crowley,  1981)  and  once  again  extremely 
efficient  algorithms  resulted  in  very  accurate 
est i mates . 

CONCLUSION 

We  have  outlined  above  several  problems  to 
which  our  spline  based  approximation  tech¬ 
niques  can  he  applied  with  great  success. 

Both  theoretical  and  numerical  findings  (some 
reported  in  the  literature  cited,  some  as  yet 
unreported  in  manuscripts)  support  our  claim 
that  these  methods  have  even  wider  applicabi¬ 
lity  than  we  have  indicated  here.  For  ex¬ 
ample,  we  are  currently  successfully  applying 
the  methods  for  estimation  of  parameters  in 
nonlinear  FDE  (Banks  and  Oaniel,  1981a)  to 
the  study  of  models  for  the  entyme  column 
reactors  as  discussed  in  (Banks,  1981), 

(Daniel,  1981).  As  one  might  anticipate  from 
the  elasticity  examples  mentioned  above,  both 
the  theoretical  soundness  and  computational 
feasibility  of  our  methods  have  been  demon¬ 
strated  for  hyperbolic  systems.  In  particular, 
we  have  successfully  developed  the  theory  and 
computational  packages  to  treat  test  problems 
in  seismic  inversion  (see  (Banks,  1981)) 
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